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3.8-Micron Photometry During the Secondary Eclipse of 
the Extrasolar Planet HD 209458b 
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ABSTRACT 

We report infrared photometry of the extrasolar planet HD 209458b during the 
time of secondary eclipse (planet passing behind the star). Observations were ac- 
quired during two secondary eclipses at the NASA Infrared Telescope Facility (IRTF) 
in September 2003. We used a circular variable filter (1.5% bandpass) centered at 
3.8 /im to isolate the predicted flux peak of the planet at this wavelength. Residual 
telluric absorption and instrument variations were removed by offsetting the telescope 
to nearby bright comparison stars at a high temporal cadence. Our results give a 
secondary eclipse depth of 0.0013 ± 0.0011, not yet sufficient precision to detect the 
eclipse, whose expected depth is ~ 0.002 — 0.003. We here elucidate the current ob- 
servational limitations to this technique, and discuss the approach needed to achieve 
detections of hot Jupiter secondary eclipses at 3.8 /im from the ground. 

Key words: infrared: techniques: photometric - eclipses - starsdndividual: HD 209458 
- planetary systems - infrared 



1 INTRODUCTION 

The passage of a transiting extrasolar planet behind its 
star during secondary eclipse has emerged as a valu- 
able technique for measuring thermal radiation emitted by 
hot Jupiters. Secondary eclipse detections of hot Jupiters 
have been accomplished using the Spitzer Space Tele scop e 
(|Charbonneau et alj|2005l : iDeming et alj|2005l f2006ah . The 
spectra of hot Jupiters are predicted to exhibit prominent 
flux peaks near 2. 2 fc 3.8 /im, where absorption by water 
vapor is minimal l|Sudarskv et all 120031 : ISeager et al. 2005; 
iFortnev et al.l 120061 ). Unfortunately, the bandpasses avail- 
able from Spitzer do not sample these flux peaks in an opti- 
mal way. Spitzer has no capabilit y at 2 /im, and th e IRAC 
filter bandpass at 3.5 and 4.5 /im (|Fazio et al.ll2004T ) overlap 
the 3.8 /im peak minimally. 

Observing hot Jupiter extrasolar planets at 2.2 and 3.8 
/im in the near future must depend on ground-based obser- 
vations. Ground-based observers have attempted both pho- 
tometry l|Snellenl 120051 : ISnellen fc Covinol 120061 ) and spec- 
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troscopy (I Richardson et all [2003) to detect the 2.2 /im 
peak using secondary eclipses. Recently, a tentative de- 
tection of OGLE-TR-1 13 at 2.2 /im was reported by 
ISnellen fc Covinol (|2006[). wit h eclipse amplitude 0.0017 ± 
0.0005. Knutson et alj |2006l ) attempted to detect the 3.8 



/xm peak (L-band) of TrES-1 spectroscopically, but no in- 
vestigations have yet reported photometry of hot Jupiters 
at 3.8 /im. At this longer wavelength, thermal background 
emission from the telescope and terrestrial atmosphere is a 
major impediment. Moreover, observations in the L band 
are significantly affected by telluric water vapor absorption, 
which is notoriously variable. In addition, photometric tran- 
sit searches at visible wavelengths have been significantly 
affected by 'red noise ' l|Pont et al.ll2006h . and this problem 
might become more severe in the infrared (IR). It is there- 
fore of interest to explore the limits of ground-based L-band 
photometry for secondary eclipse detection. 

In this paper we report L-band photometry during the 
secondary eclipse of HD 209458b. A prelim inary descrip- 
tion of these measurements was reported by IDeming et al.l 
(2006b]); here we report the final results. We do not achieve a 
detection of the planet's thermal emission, but we elucidate 
the observational limitations, and we describe improvements 
that will allow an L-band detection from the ground. 
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Figure 1. Telluric transmission near 3.8 /im, at a spectral 
resolution of X/SX = 100, based on the transmission atlas of 
iLivingston fc Wallace! l|l99lf) . shown in comparison to the stan- 
dard L' bandpass (wide dashed curve), and our CVF bandpass 
(narrow dashed curve). 

2 OBSERVATIONS 

We observed HD 209458 during two secondary eclipses on 
2003 September 9 & 16 UT, using NASA's 3-metre In- 
frared Telescope Facility (IRTF) on Mauna Kea. We also 
observed (and detected) a primary eclipse (transit) on 2003 
August 17 UT. Since t he HD 209458b tra nsit is well ob- 
served at both visible ferown et al.l l200lf ) and IR wave- 
lengths l|Richardson et al.ll2006h . we concentrate here on the 
secondary eclipse observations. However, observations on the 
transit night were very useful in analyzing the properties of 
our data (see below). 

We imaged H P 209458 using the 256x256-pixel 
NSFCam IR camera (|Shure et ail 1 19941 ). To avoid detector 
saturation from thermal background radiation, we used a 
circular-variable-filter (CVF) with a 1.5% bandpass tuned 
to 3.8 fim. This wavelength not only corresponds to a pre- 
dicted peak in the planet's thermal emission, but also has op- 
timal transmission through the terrestrial atmosphere. Fig- 
ure 1 shows the CVF bandpass in comparison to the tel- 
luric transmittance at Kitt Peak, o btained by convolving 
the at mospheric transmission atlas of lLivingston fc Wallace! 
|l99ff ) to lower spectral resolution (\/S\ = 100, compara- 
ble to CVF resolution) . Since Mauna Kea is higher and drier 
than Kitt Peak, Figure 1 is a worst case representation of 
the telluric transmission for our observations. 

It is generally not feasible to monitor comparison stars 
simultaneously with the target star in ground-based thermal 
IR photometry. Attaining the requisite wide field of view 
would require an increase in the per-pixel solid angle seen 
by the detector, to the point where the observations would 
be unavoidably saturated by thermal background. Instead, 
our observations alternated between HD 209458 and a com- 
parison star, moving the telescope between them at the most 



rapid possible cadence (17 sec for telescope motion and set- 
tling). Our imaging sequence consisted of a continuing series 
of exposures in the order T a , C a , Cb,T b ,T a , where T indi- 
cates the target star (HD 209458), and C indicates a com- 
parison star (see below). The subscripts a and b indicate 
two distinct nod positions on the detector array, separated 
by 6 arcsec (110 pixels). The nod is used for subtraction of 
thermal background (see Infrared Photometry section) . The 
observations commenced each night before the time of sec- 
ondary eclipse, and continued for as long as possible after 
eclipse. 

The ideal comparison star would have an IR bright- 
ness comparable to the target, and would lie at a close 
angular distance. Unfortunately, no ideal star exists near 
HD 209458. We must choose between fainter stars in close 
proximity, or brighter stars at greater angular distances. 
Since these are exploratory observations, we elected to try 
both approaches. On 9 Sept we monitored HD 209346 
(A2, V = 8.3, 0.2 degrees distant). To compensate for 
this star being 0.7 magnitudes fainter than HD 209458, 
we added additional comparison images to the observing 
sequence, viz.: T a ,Ca,C b ,Ca,C b ,T b ,T a , .... On 16 Sept we 
monito red HD 210483 (G1,V=7 .6, 1.7 degrees distant), fol- 
lowing [Richardson~etj2 ||2003l ). For all stars, images con- 
sisted of two co-adds of 5 sec exposures. A single 5 sec 
exposure produced about 2 x 10 s electrons of background 
radiation within our synthetic aperture for photometry (see 
below), and about 5 x 10 6 electrons due to HD 209458. 



3 INFRARED PHOTOMETRY 

Because the thermal background is intense at this wave- 
length, weighting the stellar image by the average point- 
spread-function (PSF) could in principle achieve an opti- 
mum signal-to-noise ratio (SNR) for photometry. The stellar 
PSF at this wavelength exhibits a significant component due 
to diffraction, as shown in Figure 2. Unfortunately, the PSF 
also contains a large contribution from seeing, and did not 
prove to be sufficiently stable to utilize PSF-weighted pho- 
tometry. We therefore extracted photometric intensitites for 
the star using simple aperture photometry, with a circular 
synthetic aperture of radius 20 pixels (1.1 arcsec). The best 
radius for the synthetic aperture was determined by mini- 
mizing the scatter in the photometric intensities. With radii 
significantly larger than 20 pixels, the greater background 
noise within the aperture degrades the SNR, and with sig- 
nificantly smaller radii the photometry becomes sensitive to 
seeing fluctuations. 

A significant factor in the quality of our photometry is 
telescope focus, which changes with temperature and has to 
be monitored and corrected by the IRTF observer. We found 
that both the background level and the stellar intensity var- 
ied with telescope focus. Figure 3 shows the background 
level for the HD 209458 transit observations (best night to 
illustrate this effect). Improvements in telescope focus are 
accompanied by a decrease in the measured background. 
We attribute this to less off-axis acceptance of warm radi- 
ation from the telescope structure as focus improves. The 
dependence of background on observed airmass allows us 
to deduce that ~ 80% of the background originates from 
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Figure 2. Background-subtracted image of HD 209458 at 3.8 fim, 
in reverse grayscale, showing telescope diffraction effects. The 
image core has been saturated to illustrate the weak diffraction 
wings, which are asymmetric due to optical imperfections. At this 
wavelength the PSF is a mixture of diffraction and seeing. 
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Figure 3. Thermal background intensity in our 3.8 fj,m bandpass, 
as observed on 17 August 2003. The vertical lines indicate times of 
telescope re-focusing, showing that the background level decreases 
with updates in telescope focus. 

the telescope and warm optics, the remaining ~ 20% being 
contributed by the terrestrial atmosphere. 

In principle, the stellar photometry should be indepen- 
dent of telescope focus, as long as the synthetic aperture 
is of sufficient radius to encompass slightly defocused stel- 
lar images. In practice, we find that the stellar photometric 
values change by ~ 1% when updating telescope focus. We 



found a weak correlation between intensity and the width 
of the stellar PSF (correlation coefficient ~ 0.4). Broader 
stellar images - from poorer seeing - tended to yield less 
intensity from aperture photometry. Although this correla- 
tion was dominated by seeing fluctuations, we expect that 
imperfections in focus have a similar effect. The sense of the 
correlation - less photometric intensity from broader PSFs 
- is consistent with not recovering all of the stellar photons 
within a given synthetic aperture as the image broadens. Us- 
ing extremely large synthetic apertures is impractical due to 
background noise. We therefore analyzed our data in blocks 
between focus updates, measuring HD 209458 relative to the 
comparison star. 

Although thermal IR photometry is daunting in many 
respects, there are some compensating advantages. These in- 
clude insensitivity to intra-pixel variations, and the ability 
to flat-field the detector array nearly simultaneously with 
the data acquisition. Since the thermal background will be 
spatially uniform over our 14 arcsec field of view, we use the 
thermal background for flat fielding. Within each block de- 
fined by focus updates, we compute the median image. Since 
telescope nods and pointing jitter move the star by signif- 
icantly more than its FWHM, the median image does not 
contain the star, and defines an accurate flat field calibration 
for that data block. 

We measure the intensity ratio between HD 209458 and 
the comparison star, and we do this separately for the 'a' 
and 'b' nod positions in each data block. To increase the 
precision of the comparison star photometry, we smooth the 
individual comparison measurements using a 10-point mov- 
ing average. This increases the effective time scale for our 
HD 209458 to comparison ratio to ~ 20 minutes. We spline- 
interpolate the comparison star moving average to the times 
of the HD 209458 observations, and compute the ratio. Re- 
sults from this process, for two representative data blocks, 
are shown in Figure 4. We explored the possibility of im- 
proving the precision by decorrelating the stellar intensities 
against the widths of their PSFs, but we found that the cor- 
relation was not sufficiently strong to significantly improve 
our results. 

We compute errors for each photometric measurement 
of HD 209458. Both the HD 209458 photometry and the 
comparison star photometry contribute to the total error. 
The largest source of error is background shot noise, and we 
compute its magnitude by measuring the per-pixel fluctu- 
ations in the flat-fielded background for each frame, using 
the region immediately adjacent to the photometry aper- 
ture. This noise was about two times larger than the theo- 
retical value from the square root of the electron number. 
Scaling the measured fluctuation as the square root of the 
number of pixels, we calculate the shot noise in the synthetic 
aperture due to background. Stellar photon noise and detec- 
tor read noise are negligible by comparison. We propagate 
these errors through the relative photometry, and compute 
the error in the 10-point moving average of the comparison 
stars. Comparison star and HD 209458 errors are added in 
quadrature, yielding an average per-point fractional preci- 
sion of 0.009 per single HD 209458 measurement. (In this 
paper we quote eclipse depths and errors in units of the 
stellar intensity, not in magnitudes.) This does not include 
errors due to the terrestrial atmosphere, which may not can- 
cel perfectly between HD 209458 and the comparison star. 
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Figure 4. Sample 3.8 fim photometry UT 9 Sept 2003 (upper 
panel) and UT 16 Sept (lower panel). The open diamonds show 
observations of the comparison star (HD 209346 on 9 Sept & 
HD 210483 on 16 Sept). The comparison star points have been 
normalized to 90% of HD 209458 (solid symbols), for clarity. The 
solid line is a 10-point moving average of the comparison star 
data; the dashed line scales this relation to HD 209458. Represen- 
tative error bars (±2c for clarity, and including only background 
fluctuations) are plotted over the first points in each series. Note 
the tendency for HD 209458 to track the nearby comparison star 
on 9 Sept, but not the more distant comparison star observed on 
16 Sept. 



4 RESULTS AND DISCUSSION 

Although the photometry exhibits a general decrease in in- 
tensity with increasing airmass, there are significant fluctu- 
ations on shorter time scales. Measurements of comparison 
stars are necessary for this reason. Examination of the pho- 
tometry for the nights of 9 & 16 Sept shows much better 
results using the closer, albeit fainter, comparison star (9 
Sept). This is obvious on Figure 4, where the lower panel (16 
Sept) exhibits fluctuations on time scales of ~ 20 minutes 
that are not similar between HD 209458 and the compari- 
son star. We find that this holds on 16 Sept even when fewer 
points are used in the moving average of this bright compar- 
ison star. Therefore, the data on this night are not useful for 
eclipse detection. We conclude that the relatively large an- 
gular separation between HD 209458 and HD 210483 (1.7 
degrees) is the principal factor contributing to this signifi- 
cant degree of atmospheric noise. Physically, this probably 
results from incoherence of water vapor absorption over the 
distances separating the lines of sight to these two stars. 
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Figure 5. Results for 3.8 /im photometry of HD 209458b on 
UT 9 Sept 2003. Upper panel: 331 individual observations, shown 
without error bars for clarity. Lower panel: results in bins of 0.002 
in phase, with error bars and the best-fit eclipse curve added. The 
best-fit eclipse depth is 0.0013±0.0011, not a significant detection. 



Note tha t this conclusion will no t necessarily apply to the 
results of Richardso n et al.l ([2003), since those authors used 
HD 210483 to perform spectroscopy relative to HD 209458. 

In contrast to the 16 Sept results, the data for 9 Sept 
(e.g., upper panel of Figure 4) show that fluctuations in 
HD 209458 consistently track those of HD 209346, even 
though the effective time scale of these data is ~ 20 minutes 
(10-point moving average of the comparison star). This is 
also true for the other data blocks on this night. The in- 
crease in thermal background per unit airmass was approx- 
imately two times lower on 9 Sept than 16 Sept (10% vs. 
20%). Nevertheless, we judge that the primary difference in 
our results is due to the comparison star selection: a nearby, 
albeit faint, comparison star is more useful than a brighter 
but more distant one for photometry in the L-band. Cer- 
tainly this is true in conditions that are less than optimum. 

Figure 5 (upper panel) plots the 331 individual obser- 
vations from 9 Sept versus the plan e t's or bital phase, using 
the ephemeris from iRnutson et al.l (|2007l ). These data are 
normalized to an average of unity. The lower panel aver- 
ages the data in bins of width 0.002 in phase, and adds 
error bars calculated as discussed above. The best-fit eclipse 
depth (0.0013 ± 0.0011) does not have sufficient precision 
for detection, but our re sults are simila r to a r esult recently 
reported for TrES-1 by iKnutson et al] l|2006l ). who used a 
spectroscopic technique at 3.8 /xm. 

The reduced chi-squared of the fit in the lower panel 
of Figure 5 is 1.96, consistent with a per-point scatter of 
0.014 (upper panel), about 50% larger than the errors from 
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background shot noise. Additional error can arise from at- 
mospheric noise that is not common to both stars on the 
spatial and temporal scales of these observations. This 'red 
noise' (amplitude ~ 0.01) can be reduced by using a broader 
optical bandwidth. New IR cameras such as NSFCam-2 at 
the IRTF can utilize the full L' band without saturating. 
This will increase the SNR by a factor of ~ 3.5 per image 
(SNR proportional to square root of optical bandwidth). 
Having greater SNR per frame should increase the preci- 
sion of fainter comparison stars; for example, HD 209346 
should improve from 0.021 to 0.006 precision in 10 seconds 
of integration. A more significant effect of a broader optical 
bandwidth is that it will decrease the temporal bandwidth 
of the red noise, by an order of magnitude. This follows be- 
cause the time to define the target-to-comparison-star ratio 
to a given precision shortens as the square of the SNR per 
image, and we would not need to utilize a 10-point moving 
average for the comparison star. 

Atmospheric noise power often increases as the inverse 
of the sampling frequency ('1/f noise). We have calculated 
the improvement expected in both the background shot 
noise and the atmospheric 1/f noise using the full L' opti- 
cal bandwidth. This calculation assumes that any high fre- 
quency contribution from the terrestrial atmosphere - sim- 
ilar to scintillation - does not dominate on ~ 10 to 30 sec- 
ond time scales at this wavelength. The result indicates that 
our 9 Sept final errors will decrease by a factor of two. A 
secondary eclipse of depth 0.002 in HD 209458 would be 
detected to 3.6<r significance. Moreover, the greater planet- 
to-st ar contrast of the r ecently discovered HD 189733 sys- 
tem dBouchv et al.ll2005l ) , coupled with its brighter apparent 
magnitude, will lead to even more favorable detectability. 
We conclude that ground-based detection of photons from 
extrasolar planets is possible at 3.8 /im. 
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